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Abstract 

Fusion data for ^^C+^^C, ^^C+^^C and ^^C+^^C are analyzed by coupled-channels calculations 
that are based on the MSY+repulsion, double-folding potential. The fusion is determined by 
ingoing-wave-boundary conditions (IWBC) that are imposed at the minimum of the pocket in the 
entrance channel potential. Quadrupole and octupole transitions to low-lying states in projectile 
and target are included in the calculations, as well as mutual excitations of these states. The 
effect of one-neutron transfer is also considered but the effect is small in the measured energy 
regime. It is shown that mutual excitations to high-lying states play a very important role in 
developing a comprehensive and consistent description of the measurements. Thus the shapes of 
the calculated cross sections for ^^C+^^C and ^^C-|-^^C are in good agreement with the data. The 
fusion cross sections for ^■^C+^^C determined by the IWBC are generally larger than the measured 
cross sections but they are consistent with the maxima of some of the observed peak cross sections. 
They are therefore expected to provide an upper limit for the extrapolation into the low-energy 
regime of interest to astrophysics. 

PACS numbers: 24.10.Eq, 25.60.Pj, 26.30.-k, 26.50.+X 
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I. INTRODUCTION 



The fusion of carbon nuclei is an important reaction in the description of type la su- 
pernovae and other astronomical events in the cosmos like the superburst of an accreting 
neutron star. There is, unfortunately, a very large uncertainty in the predicted fusion cross 
sections for ^^C+^^C that are needed in a stellar environment partly because the cross 
sections are difficult to measure down to the low energies that are of interest, and partly 
because the data contain strong resonance structures that make it difficult to extrapolate the 
measured cross sections to low energy. In order to get some constraints on the extrapolation 
it is instructive to analyze the existing fusion data for ^'^C+^^C by Trentalange et al. l2| and 

n n f] 

Charvet et al. [3|, and also of the ^^C+^'^'^C data by Dayras et al. |^] and Notani et al. [5|], 
because these data do not exhibit strong resonance features. 

The ^^C+^'^C fusion data [2l] are analyzed by coupled-channels calculations that include 
couplings to the one-phonon quadrupole and octupole excitations in projectile and target, 
as well as mutual excitations of these states. The effect of one-neutron transfer, which can 
proceed to the ^^C+^^C mass partition with positive Q value is also studied, and so is the 
effect of the neutron exchange with zero Q value in ^^C+^'^C reactions. 

The fusion cross sections are determined by ingoing-wave-boundary conditions (IWBC) 
that are imposed at the minimum of the pocket in the entrance channel potential. The 
calculated cross sections defined this way are fairly smooth as functions of the center-of- 
mass energy and they are well suited to analyze the fusion data for ^^C+^'^C and 
i2q_|_i3q 2^]. The ^^C-|-^^C fusion data, on the other hand, contain a lot of structures or 
resonances, and it is beyond the scope of this investigation to try to reproduce these data 
in detail. It is, however, of great interest to see how the calculated cross sections, obtained 
from the IWBC, compare to the data and how they possibly can be used to put constraints 
on the extrapolation to very low energies. 

The coupled equations are solved using either a standard Woods-Saxon \^ or the 
M3Y-I- repulsion, double-folding potentials [tJ. The coupled-channels effects on the calcu- 
lated fusion cross sections are relatively modest compared to the large enhancement of 
several orders of magnitude that are commonly seen in calculations of heavy-ion fusion reac- 
tions. For the carbon systems, the enhancement of fusion compared to the no-coupling limit 
is typically a factor of 2 at energies far below the Coulomb barrier. However, it is a challenge 
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to explain the data in detail because the coupled- channels calculations are sensitive to chan- 
nels that have rather high excitation energies. One problem is that these channels are closed 
at low beam energies but that problem can be solved by imposing the correct, decaying state 
boundary conditions at large distances between the reacting nuclei jsj. Another problem is 
that the nuclear structure of high-lying states is sometimes very uncertain. 

The M3Y+repulsion potential was introduced in Ref. I?] to explain the hindrance that 
has been observed in the fusion of many heavy- ion systems at very low energies 9|. The 
hindrance phenomenon was first observed as a suppression at very low energies of ^^Ni+^^Y 
usion data compared to calculations that were based on a standard Woods-Saxon potential 
lo| . A simple explanation of the phenomenon is the existence of a shallow pocket in the 
entrance channel potential which forces the fusion cross section to vanish as the center- 
of-mass energy approaches the minimum energy of the pocket [3]. Important issues have 
been whether the fusion hindrance also occurs in light-ion systems, and how it will affect the 
extrapolation of measured cross section to the low energies that are of interest to astrophysics 
11 1 . The evidence for quasi-molecular resonances, observed in the elastic scattering |l2| and 



fusion reactions 
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14| of ^^C+^^C are usually explained as resonances in a shallow two- 
body potential between the reacting nuclei. We will show that the analysis of the ^^C+^^C 
and ^^C+^'^C fusion data supports the idea of a fusion hindrance and the existence of a 
shallow pocket in the entrance channel potential. 

The fusion cross sections reported in Refs. 2|-|^ are based on measurements of the char- 
acteristic 7-rays emitted from some of the evaporation residues that are produced, mostly 
those associated with the proton, neutron, and a decay. The total fusion cross sections were 
obtained with the aid of statistical model calculations which is a major source of the sys- 
tematic uncertainty. The systematic error is difficult to estimate but it can be quite large. 
For example, the systematic error of the absolute cross sections quoted in Ref. % is ±15%. 
Since some of the systematic error concerns the overall normalization of the measured cross 
sections, and not so much the shape (i. e., the energy dependence of the cross section), it is 
of interest to adopt an adjustable overall normalization when analysing the data. 

The ingredients of the coupled-channels calculations are presented in the next section. 
The analysis of the ^^C+^^C fusion data is discussed in section III, where the optimum 
parameters of the M3Y+repulsion potential, such as the radius parameter of ^^C and the 
diffuseness associated with the repulsive term, are determined. These parameters are used 



3 



together with a parameterization of the experimental density of ^^C as input to the coupled- 
channels calculations of the fusion cross sections for ^^C+^^C and ^^C+^^C which are pre- 
sented in sections IV and V, respectively. The systematics of the low-energy fusion cross 
sections for the 3 systems of carbon isotopes is discussed in section VI. Finally, the conclu- 
sions of this work are presented in section VII. 



II. DETAILS OF THE CALCULATIONS 



The ingredients in the coupled- channels calculations include the ion-ion potential, the 
nuclear structure input, the strength and Q- value of the transfer reactions, and the definition 
of the fusion cross section. The coupled-channels technique that is used has been applied 
previously and is described, for example, in Refs. j3, M, H, All of the details will 
therefore not be repeated here. It is emphasized that the calculations are performed in the 
rotating frame approximation 



16] which makes it feasible to include the most important 
reaction channels in the calculations. In this approximation, the magnetic quantum number 
M of the entrance channel is conserved. In reactions between nuclei with 0"*" ground states, 
this implies that the M quantum number remains zero. In reactions with odd nuclei one 
would have to repeat the calculations for each initial value of M and the average fusion cross 
section should be compared to the measurements. In the fusion of ^^C+^^C, for example, 
the ground state of ^^C is a 0^ state and ^^C has a l/2~ ground state. However, the cross 
sections for M = -1/2 and +1/2 are identical so it is sufficient to calculate the cross section 
once. 



A. Densities and ion-ion potentials 

Besides the commonly used Woods-Saxon potential, the ion-ion potential that will be used 
in this paper is the so-called M3Y+repulsion potential It consists of the conventional 
M3Y potential and a repulsive term which is discussed below. The M3Y double-folding 
potential is defined as 

[/„(r) = J dri dr2 pi(ri) p2(r2)^^M3y(r + ra - ri), (1) 

where VM3Y{f) is the M3Y effective nucleon-nucleon interaction derived from the Reid po- 
tential 18|. The densities pi{r) of the ^^C and ^^C carbon isotopes are parametrized in terms 



of the symmetrized fermi functions defined in the appendix of Ref . [17| . One advantage of 
this parametrization is that the mean square radius is given by the simple expression, 



in terms of the radius R and the diffuseness a. 

The diffuseness of the densities for both nuclei is set equal to 0.44 fm so the only param- 
eters that need to be specified are the matter radii of the reacting nuclei. This choice of the 
diffuseness is in good agreement with the tail of the measured point-proton density of ^^C 
{19I which is illustrated in Fig. [H The calculated point-proton distribution has the radius 
R = 2.155 fm consistent with the measured rms charge radius. The radii that reproduce the 
point-proton (pp) rms radii of ^^C and ^'^C, extracted from the measured root-mean-square 
(rms) charge radii [2^, are shown in the last few lines of Table I. The quoted radius for ^^C 
is expected to be a realistic estimate of the matter radius of ^^C and it will therefore be used 
in calculations of the M3Y double-folding potentials. 

The matter radius of ^'^C is possibly larger than the radius of ^^C because of the valence 
neutron. That is what one would expect by adding a neutron to an inert ^^C core nucleus. 
A simple estimate based on a single-particle, Woods-Saxon potential gives a mean square 
radius of (r^) = 8.655 fm^ for the pi/2 valence neutron. Combined with the point-proton 
rms radius of ^^C, {r'^)i2 = 5.462(9) fm^, one can now estimate the mean square matter 
radius (r^)i3 of ^^C by 



The result is a rms matter radius of 2.378 fm. The associated radius of a fermi-function 
distribution with diffuseness a = 0.44 fm is 2.228 fm according to the last line of Table I. 
This is a useful reference value for the discussion below. 

A critical part of the M3Y+repulsion, double- folding potential is the repulsive part 
which is generated by a contact interaction, f,.5(r). The densities that are used in calculating 
the repulsive, double-folding potential have the same radius as those that are used in the 
calculation of the M3Y double-folding potential but the diffuseness is chosen differently 
(usually smaller.) The strength Vr of the repulsion, on the other hand, is adjusted for a given 
value of ttr so that the nuclear incompressibility K = 234 MeV is produced (see Ref. [7| for 
details.) With this constraint there are three parameters that must be specified before the 



{r') = liR' + l{7ran 



(2) 




(3) 
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M3Y+repulsion potential can be calculated, namely, the matter radii of the two reacting 
nuclei, and the diffuseness parameter associated with the repulsion. 

The parameters that give the best fit to Trentalange's ^'^C+^'^C fusion data are the ^^C 
radius R = 2.28 fm and the diffuseness = 0.31 fm. How these parameters were determined 
is described in detail in section III.B. The radius R = 2.28 fm is slightly larger than the 2.228 
'm that was estimated from Eq. ([3]). That may not be unreasonable as discussed in Ref. 



2l| because using a larger radius in a coupled-channels calculation is a way to compensate 



for the dynamic polarization of states that are not included explicitly in the calculation. 
The discrepancy between the estimated rms radius, Eq. ([3]), and the value extracted from 
the fit to the fusion data is only 1.2 %. 

The entrance channel potential determined by the best fit parameters is shown by the 
solid curve in Fig. |2JA.. Also shown in Fig. [2] are the potentials that will be applied in 
calculating the fusion of ^^C+^^C and ^^C+^^C. They were obtained by setting the radius 
of ^^C equal to the radius R = 2.155 fm of the point-proton distribution of ^^C (see Table 
I), and the diffuseness associated with the repulsion was set to = 0.31 fm. 

The pure M3Y double-folding potentials shown in Fig. [2] are unphysical because they 
are much deeper for overlapping nuclei than the ground state energy of the compound 
nucleus. The (blue) dashed curve in Fig. [2JA. is the entrance channel potential for a standard 
Woods-Saxon (WS) potential with the depth Vq = -39.31 MeV, radius R^s = 5.268 fm, and 
diffuseness a = 0.63 fm. Here the radius was adjusted to optimize the fit to Trentalange's 
fusion data as discussed in section HI. A. This potential is deeper than the M3Y-|-repulsion 
potential and has a minimum that is close to the energy of the compound nucleus ^^Mg. 

The three entrance channel potentials shown in Fig. |2JA., which are based on the M3Y, 
the M3Y-|-repulsion and the Woods-Saxon potential, produce essentially the same Coulomb 
barrier height. Vc^ ~ 6 MeV. The most important features of the M3Y-|-repulsion potential 
are that it produces a shallow pocket and a relatively thick barrier. These features were 
utilized in Ref. [tI to explain the hindrance of fusion that has been observed in the fusion 
of many medium mass systems at extreme subbarrier energies. The analysis of the ^^C-|-^^C 
fusion data discussed in the next section confirms that a shallow pocket in the entrance 
channel is a general feature that helps explain the energy dependence of fusion data. 
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B. Nuclear couplings 



The nuclear couplings that excite the low-lying states of projectile and target are derived 
from the macroscopic description of surface excitations j^. In this description, the surface 
of nucleus i is parametrized as 

i?.(f) = i?f)(i + E««y,;(f)), (4) 

where are the (static or dynamic) deformation amplitudes of the nucleus, and r specifies 
a spatial direction. The distortion of the nuclear surface, 

is an operator that can excite and de-excite the nucleus (see chapter II. 4 of Ref. Q].) The 
excitation can be vibrational, in which case the intrinsic hamiltonian is a harmonic oscillator, 
or it can be rotational, in which case axf^ = P\ D^Q(e), where Px is the static deformation 
parameter and e is the direction of the symmetry axis of the deformed nucleus. 

In the macroscopic description of heavy- ion reactions gI, the nuclear potential between 
two colliding nuclei is parametrized in terms of the surface distortions 6Ri as follows 

V^{r - - - 6R,{r) - 5i?2(-r)), (6) 

where f is the direction of the center of mass distance between projectile and target. In the 
rotating frame approximation, which is a simplifications that is commonly used and which 
will be used in coupled-channels calculations, the direction of f defines the z-axis, so that 



^fl,(f)=<i:4i> J!^. (7) 



In the model of heavy-ion fusion reactions described in Ref. 16|, the nuclear potential is 
expanded up to second order in the surface distortions, 

f/^(r, 5Ru 6R2) = Vn- ^{SRi + 5R2) 

ar 



+ 



1 d^VN 



2 dr^ L^5/?i + (5i?2) -(0|((5i?i + 5i?2)10)J, (8) 

where the argument of Vat and its derivatives is r — R^'^ — -R2°'*- The second-order term in 
this expression has been renormalized so that the ground state expectation value of that 



term is zero. The ground state expectation value of the first-order term (proportional to 
6R1+6R2) will also vanish if the ground states of the two reacting nuclei are O"*" states. The 
ground state expectation of the entire expression, Eq. (jS]), is therefore identical to the 'bare' 
interaction, i. e., 

(0|f/^(r,5i?i,5i?2)|0) = V^ir - - 4°^). (9) 

This implies that the M3Y+repulsion potential which is calculated using the ground state 
densities of the reacting nuclei can be used as the bare interaction, Vn- 

The diagonal matrix element of the first-order term in Eq. ([8]) can be non-zero, for 
example, in the 2+ state of a deformed nucleus. The calculation of diagonal and off-diagonal 
matrix elements is outlined in the appendix. Details of how to calculate the matrix elements 
of the second-order couplings in Eq. ([8]), both for vibrational and rotational excitations, are 



given in Ref. 16 1. 



C. Nuclear structure input 

The nuclear structure input to the coupled-channels calculations is listed in Table II. 
The nuclei ^^C and ^^C are both considered deformed with oblate quadrupole shapes. The 
connection between the deformation parameter Px and the measured B{E2) values is sum- 
marized in the appendix. For ^^C one can extract the quadrupole deformation parameter P2 
= 0.57(2) from the measured strength of the Of — t- 2+ transition. The value of this deforma- 
tion parameter produces an intrinsic quadrupole moment of Qq = -19.5 fm^ (assuming an 
oblate shape). The quadrupole moment of the 2+ state is therefore Q2 = -2/7Qo = 5.57 fm^. 



23j. Also shown in the Table is the 



which is consistent with the measured value of 6 ± 3 fm^ 
structure information about the — transition which is included in the calculations as 
part of the two-phonon quadrupole excitation of ^^C. 

Values of quadrupole deformation parameter (32 of ^'^C, extracted from Eq. ( |2Ti) of the 
appendix and the known strengths of the quadrupole transition from the l/2~ ground to 
the 3/2~ and the 5/2~ excited states, are shown in the last column of Table II. They are 
seen to be almost identical, which indicates that ^^C is a fairly good rotor. Assuming the 
deformation parameter /32 = 0.495 and an oblate shape, one obtains the intrinsic quadrupole 
moment Qo = -17.9 fm^. 

The strength of the octupole transition in ^^C quoted in Table II appears to be very large 
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and almost as large as the octupole strength in ^^C. This is misleading because what matters 
is the off-diagonal matrix element of the octupole amplitude, which is given by Eq. f[T9l) of 
the appendix. In ^^C the matrix element is 

(300|a3o|000) = (10) 

whereas in ^^C it is (with M = K= 1/2) 

.5 11, ,111, . /2 11 51 2 /3 /33 .... 

The matrix element of the octupole amplitude is therefore reduced in ^"^C by the factor 



3/7 ~ 0.655 relative to the expression Eq. flTOj) for ^^C. 

The coupled-channels calculations that are performed are similar to those presented in 



Ref. 



16| for the fusion of Al with different germanium isotopes. In addition to the nuclear 



couplin gs d erived from the nuclear interaction, Eq. ([S]), the Coulomb interaction (Eq. (2) 



of Ref. 



16| ) is included in the calculations to first order in the deformation amplitudes. 



D. Transfer reactions 

The one- neutron transfer reactions that will be considered all involve the pi/2 orbit, both 
in the initial and final states. One example is the ground state to ground state transfer 
reaction, ^^C(^'^C,^^C)^'^C, which has zero Q value. The spectroscopic factors for the initial 
and final state l/2~ state in ^'^C were set equal to 0.75(10), which is the value recommended 
in Table III of Ref. |24|. 

The other example is the one-neutron transfer from ^'^C+^^C to the ground states of 
i2q_|_i4q ^Jiic]^ j^as a Q-value of +3.2 MeV. The spectroscopic factor for the initial state is 
the same as above. It is set equal to 2*0.75 in the calculations because the transfer can take 
place from either of the two ^^C nuclei in the entrance channel. The spectroscopic factor for 



the pi/2 orbit in the 0+ ground state of ^"^C is 1.63(33) according to table III of Ref. [24 1. 



The transfer form factors that are used are the so-called Quesada single-particle form 



factors [25] which in the past turned out to be fairly realistic |l7l . \2&\ . They have here 
been calibrated against calculations performed with the computer code Ptolemy {2^ of the 
one-neutron transfer cross in ^^C+^^C collisions at a 5.2 MeV center-of-mass energy. 
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E. Calculation of fusion cross sections 



The fusion cross section is calculated from the ingoing flux which is determined by the 
ingoing-wave-boundary conditions (IWBC) that are imposed at the minimum of the pocket 
in the entrance channel potential. This approximation ignores the internal structure of the 
combined di-nuclear or compound system and the calculated cross sections are usually fairly 
smooth functions of energies, at least at energies near and below the Coulomb barrier. 

At energies far above the Coulomb barrier, there are sometimes numerical problems in 
the coupled-channels calculations which can cause an erratic behavior of the calculated 
cross section as a function of energy. One can overcome these problems by applying a weak 
imaginary potential Q 



which acts near the position R„ of the pocket in the entrance channel potential. The fusion 
cross section is then defined as the sum of the absorption cross section and the cross section 
for the ingoing flux. A weak imaginary potential will be applied in the following when the 
calculated fusion cross section is shown in a linear plot at energies far above the Coulomb 
barrier. The parameters that are used are Wq = -5 MeV and Ow = 0.5 fm. No imaginary 
potential will be used in the calculations at energies below the Coulomb barrier. 

The 1/2" ground state spin of ^^C causes some concern when calculating the scattering 
and fusion of ^'^C+^'^C. The fermionic nature of ^^C requires that the total wave function 
for ^'^C+^'^C be anti-symmetric. By coupling the intrinsic l/2~ ground state spins of the 
reacting nuclei one obtains a total spin of S=0 (anti-symmetric) and S=l (symmetric). The 
wave function for the relative motion of the two nuclei must therefore be symmetric for S=0, 
i. e., it consists of even partial waves, and it must be antisymmetric for S=l, consisting of 
odd partial waves. The fusion cross section is therefore calculated as the weighted sum 



It turns out that the contributions to the fusion cross section from odd and even L- 
values are almost identical. The biased weighting of the two contributions in Eq. ( !T3|) does 
therefore not have much effect. In calculations of the fusion and elastic scattering of two 
identical spin zero ^^C nuclei, one should only consider the contribution from even partial 
waves which is doubled because of the symmetry of two identical particles. In reactions of 




(12) 




(13) 



10 



i2q_|_i3q^ on the other hand, the contribution from even and odd partial waves have the 
same weight. 

III. ANALYSIS OF THE ^'^C+^^C FUSION DATA 

In this section the ^'^C+^^C fusion data by Trentalange et al. 2I and Charvet et al. 
are analyzed in terms of coupled-channels calculations. The channels that are included are 
enumerated in Fig. 121 The (red) solid lines represent the elastic channel and the channels 
associated with the excitation of the states shown in Table II, in both projectile or target. 
The calculation based on these channels is referred to as the one-phonon (phi) calculation 
and has a total of 7 channels. 

The (blue) dashed lines in Fig. [3] are the 6 mutual excitations of the lowest quadrupole 
and octupole excitations in projectile and target. These include the four mutual states 
(/i,J2), where Jj = 3/2^ or 5/2+, i=l,2, and Ji and I2 belong to different nuclei, and 
the two (3/2", 5/2+) mutual excitations that belong to the same nucleus. Together with 
the seven channels of the one-phonon calculation that gives a total of 13 channels. This 
calculation is referred to as the mutual excitation calculation or Chl3. 

The calculation Chl3 will be compared to the one-phonon calculation (phi) described 
above and the no-coupling limit which has only one channel. There are other excitations, 
for example, the two-phonon excitations which are poorly known and mutual excitations 
that involve the 5/2" state but they are all ignored in the following. 

It can be very difficult to see small discrepancies between measured and calculated cross 
sections, in particular, when they are plotted on a conventional logarithmic scale. It is 
therefore useful to use other representations that amplify certain features of the comparison, 
such as the S factor for fusion, which emphasizes the behavior at low energies. Another 
representation, which is very useful when coupled-channels effects are modest, is the so-called 
enhancement factor, which was used in the early days when the enhancement of subbarrier 
fusion was first discovered. It is defined as the ratio of a cross section relative to the cross 
section calculated in the no-coupling limit, and it will be used both for measurements as well 
as coupled-channels calculations, as long as it is indicated which no-coupling limit is used as 
a reference. Finally, the comparison of data and coupled-channels calculations will also be 
made in terms of the ratio of the measured and calculated cross sections. If the measured 
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and calculated cross sections agree, the ratio should be close to one. 

The one-neutron transfer discussed earlier is included as an independent degree of free- 
dom. That implies that the transfer can take place from any excited state and it is calculated 
with the same form factor which describes the ground state to ground state transition. The 
full calculation (with 13 excitation channels mentioned above and the one-neutron transfer) 
will therefore contain 26 channels and is referred to as the Ch26 calculation. 

The analysis of the data is performed in terms of the per data point, x^/^? which is 
minimized by adjusting certain parameters of the ion-ion potential. The overall normaliza- 
tion of the calculated cross sections is also adjusted by a multiplicative factor Sc, in order to 
improve the fit to the data. The motivation for this adjustment is the large uncertainty in 
the absolute normalization of the measured cross sections that was mentioned in the intro- 
duction. The shape (i.e., the energy-dependence) of the measured cross section is assumed 
to be more accurately determined. The uncertainty in the data analysis will therefore in- 
clude the statistical uncertainty and an adopted systematic error of only 3%, which is much 
smaller than the 15% experimental error 2]- 

A. Application of Woods-Saxon potential 

The first set of calculations are based on a standard Woods-Saxon potential with typical 

n 

parameters from Ref. p. The diffuseness of the potential is set to a = 0.63 fm and the 
depth is -39.31 MeV. The radius of the potential, = 5.268 fm, was adjusted to optimize 
the fit to the data of Ref. 2| . A comparison to the data is made in Fig. H] in terms of the 
enhancement of the cross sections with respect to the no-coupling limit. The enhancement 
of the measured cross sections is seen to deviate from the enhancement of the coupled- 
channels calculation (Ch26) by up to 20%. Moreover, the effects of mutual excitations and 
one-neutron transfer are seen to be modest by comparing to the no-transfer calculations 
(Chl3) and the one-phonon calculations (phi). It is very unlikely that the agreement with 
the data can be improved much by adjusting the strengths of the couplings to the excitation 
or transfer channels. It is shown in the next subsection that a much better agreement with 
the data is achieved by applying the M3Y-|-repulsion instead of the Woods-Saxon potential. 

The agreement with the data in Fig. H] cannot be improved much by scaling the calculated 
cross sections with an adjustable scaling factor Sc in the data analysis. The best fit for Sc = 
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1 has a "x^ jN = 8.7. It can be reduced to x^/^ = 6.8 by adjusting both the scahng factor 
Sc and the radius Rws of the Woods-Saxon well. The values that give the best fit are Sc = 
0.885 and R^s = 5.345 fm. 

Another observation in Fig. H] is that the data sets of Trentalange et al. ^ and Charvet 
et al. [3j] differ by up to 10% in the overlapping energy regime. However, this is not a serious 
problem because the systematic error of both experiments is about 15%. 



B. Application of the MSY+repulsion potential 

One advantage of the M3Y+repulsion potential is that it has an additional adjustable 
parameter besides the radius R of the density, namely, the diffuseness parameter of the 
density that is used in calculating the repulsive part of the potential. This parameter controls 
the depth of the pocket and the thickness of the barrier in the entrance channel potential. 

The determination of the best fit parameters is illustrated in Fig. [5l For a fixed value of 
Qr, the fusion cross sections were calculated for different values of the radius R of ^'^C. In 
each case the scaling factor Sc was adjusted to give the best jN . The results are shown by 
the dashed curves. The minimum /N for each dashed curve defines the solid curve, and 
the minimum of this curve defines the absolute best fit (ABF) to the data. The parameters 
of the minimum are = 0.31 fm, the matter radius R = 2.28 fm, and the scaling factor Sc 
= 0.843. They are compared in Table 1 to other parameters that were obtained for fixed 
choices of the radius of ^^C. 

The x^ /N of the best fit to Trentalange's data is small according to Table 1, x^ /N =1.0, 
although the adopted systematic error of the analysis was set to only 3%. The scaling 
factor of the best fit, Sc = 0.843, is not unreasonable but it is slightly extreme because the 
systematic uncertainty in the absolute cross section of the measurement was estimated to 
be 15% 2I. It is therefore of interest to discuss whether the other parameters of the best fit 
are realistic. The diffuseness ar associated with the repulsion, ar = 0.31 fm, falls below the 
range of values ar = 0.4 - 0.43 fm that have been obtained in the analysis of fusion data for 



other symmetric heavy-ion systems 



2l|. 



The ^^C radius of the best fit, R = 2.28 fm, produces an rms matter radius of 2.407 fm 
according to Table 1 which is slightly larger than the rms matter radius of 2.378 fm that 
was estimated in Eq. This result is not unreasonable as discussed in connection with 
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an analysis of the fusion data for ^^Ca+^^Ca [21]. There it was pointed out that a larger 
radius may simulate the effect of the dynamic polarization of excited states which are not 
included explicitly in the coupled-channels calculations. 

If the scaling parameter Sc is kept fixed at Sc=l in the data analysis, the best fit is 
achieved with a smaller ^^C radius, namely, R = 2.17 fm. This value is almost as small as 
the radius for the point-proton density distribution of ^'^C. The fit is poor, with a /N = 
2.75 as shown in the first line of Table I. The value of the diffuseness parameter ar which 
produces the best fit in this case is = 0.33. 

The ratios of the measured and calculated cross sections discussed above are illustrated 
in Fig. O The ratio of Trentalange's data I2I to the best fit (BF) solution with radius R 
= 2.17 and = 1 is shown by the open symbols. The ratio of the same data set to the 
absolute best solution (ABF) with radius R = 2.28 fm is shown by the solid (red) points. 

n 

Also shown at higher energies is the ratio of Charvet's data [3] to the same solution (with 
radius R = 2.28 fm.) It is seen that the ratios for the absolute best solution are essentially 
constants (~ 0.84 for Trentalange's data, and ~ 0.92 for Chatvet's data up to 10 MeV). The 
two data sets differ by up to 10% in the overlapping energy regime. However, this is not 
a serious problem as mentioned earlier because the systematic error of both experiments is 
about 15%. 

The ratio for best fit (BF) solution for Sc=l exhibits some energy dependence in Fig. [61 
It is of course possible that the data do contain some structures (similar to those observed in 
the ^^C+^^C fusion data) that should not be accounted for by the type of coupled-channels 
calculations performed here. The absolute best fit (ABF) solution, obtained with the ^^C 
radius R = 2.28 fm and Sc = 0.843, could therefore be misguided or misleading. On the 
other hand, the absolute best fit does have some attractive features as argued above: the 
quality of the fit is excellent, and the extracted rms matter radius is only 1.2% larger than 
expected. The absolute best fit requires a scaling of the calculation, Sc = 0.843, which is not 
unreasonable considering the large 15% systematic error of the experiment. The parameters 
of the ABF are therefore adopted in the following. 
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C. Results of the analysis 



In order to illustrate what makes it possible to reproduce the energy dependence of 
Trentalange's data so well when using the M3Y+repulsion potential, we show in Fig. [7] the 
enhancement of the calculations and the data relative to the no-coupling calculation. In 
contrast to the results for the Woods-Saxon potential shown in Fig. HI there is now a strong 
sensitivity to mutual excitations. This is caused by a larger second derivative of the ion-ion 
potential (stemming from the onset of the strong repulsion), and consequently, a stronger 
quadratic coupling of the entrance channel to mutual excitations. 

An interesting observation in Fig. [7] is that the neutron transfer does not have much 
influence on the calculated cross section in the energy regime of the experiment. This is 
seen by comparing the calculations labeled Chl3 and Ch26. The influence shows up at 
energies below 3 MeV, so it will have an impact on the extrapolation of the cross sections 
to lower energies. Another interesting feature in Fig. [7] is the structure of the data at low 
energies. The small peaks at 3.6 and 4.2 MeV, for example, could be remnants of similar 
structures observed in the fusion data for ^^C+^^C. 

The calculated S factors for the fusion of ^'^C+^^C are compared in Fig. [H]to the data of 
Ref. {21. The data have here been divided by the scaling factor Sc = 0.843 which optimizes 
the agreement with the full calculation labeled Ch26. The calculation that is based on 
one-phonon excitations only (phi) does no reproduce the energy dependence of the data 
between 3 and 6 MeV. The calculation that includes mutual excitations but no transfer (the 
Chl3 calculation) is seen to reproduce the corrected data very well. Finally, the additional 
coupling to the one-neutron transfer produces some enhancement in the full Ch26 calculation 
but that occurs at energies below the range of the experiment. 

The calculated cross sections obtained with all 26 coupled channels and the parameters 
determined in the previous subsection are compared in Fig. [9] to the data of Ref. j^. 
The calculation exceeds the data on average by only 4%, which is small compared to the 
15% systematic error of the experiment. The coupled- channels effects are fairly modest at 
high energies; this can be seen by comparing to the no-coupling calculation which is shown 
by the thick dashed curve. The results of a maximum angular momentum cutoff in the 
coupled-channels calculations are shown for Lmax = 2-12. 
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IV. PREDICTIONS OF THE FUSION OF i^C+^^C 

Having determined the radius of ^^C and the diffuseness a^. associated with the repulsive 
part of the double-folding potential in the previous section, one can now predict the ion-ion 
potential for ^^C+^^C, provided the density of ^^C is known. Here the parameters of the 
point-proton density of ^^C given in Table 1 are adopted. The predicted M3Y+repulsion 
potential is shown in Fig. |2j3- It is similar to ^^C+^^C potential shown in Fig. |2j^. The 
main difference between the two figures is that the ground state of the compound nucleus 
^^Mg is located at a higher energy than ^^Mg because of the unpaired neutron in ^^Mg. 

The excitation energies that are considered in the calculations of the fusion cross sections 
for ^^C+^'^C are shown in Fig. [101 The solid (red) lines represent the excitations shown 
in Table I which together with the elastic channel give a total of 7 channels. One of these 
states is the 0^ Hoyle state of ^^C which is included because its coupling to the 2"*" is known. 

The (blue) dashed lines in Fig. [10] are the 6 mutual excitations that are considered. They 
include the (2+,3~) mutual excitation in ^^C, the (3/2~, 5/2+) mutual excitation in ^^C, 
and the four (/i,/2) mutual excitations in projectile and target with Ji = 2+ or 3^ and 
I2 = 3/2^ or 5/2+. Two-phonon excitations (except the 0^ Hoyle state which has already 
been mentioned) and mutual excitations that involve the 5/2" state in ^^C are ignored. The 
total number of excitation channels is therefore 7+6 = 13. When the one-neutron transfer 
with zero Q value is included in the calculations as described in Sect. II. D, the number of 
channels is doubled. The full calculation, which includes the one-neutron transfer and all of 
the 13 excitation channels shown in Fig. [TD]will therefore have 26 channels and it is referred 
to as Ch26. 

A major problem, which causes some uncertainty in the predictions made by the coupled- 
channels calculations, is the high excitation energy and the large strength of the octupole 
transition in ^^C (see Table 1.) A large excitation energy implies decaying boundary con- 
ditions at large separations of projectile and target instead of the scattering boundary con- 
ditions, and this can cause numerical problems. A large octupole strength implies a strong 
sensitivity to mutual excitations that involve the 3~ state and some of these excitations are 
poorly known as explained below. The octupole excitation of ^'^C, on the other hand, is not 
a problem because the excitation energy of the 5/2+ state is much lower than in ^^C, and 
the effective octupole strength is reduced as shown in Eq. (fTTl) . 
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The enhancement of the ^^C+^^C fusion relative to the no-couphng hmit is shown in Fig. 
[11] for different coupled-channels calculations. The full calculation Ch26 has a very strong 
peak near 4.7 MeV. The calculation Ch24, which excludes the mutual (2+,3^) excitation 
in ^^C, has a much weaker peak. There are two reasons why the mutual (2+, 3~) excitation 
in ^^C has such a large influence on subbarrier fusion. One reason is the large octupole 
strength of ^^C. The other reason is that the direct coupling of the ground state to mutual 
excitations is governed by the second derivative of the ion-ion potential, according to Eq. 
dH]), and this quantity is particularly large for the M3Y-I- repulsion potential as pointed out 
in section III.C. 

It is very interesting that the enhancement of the data % |5| with respect to the no- 
coupling limit also exhibits a peak structure in Fig. [TTJ The peak is located near 4.3 MeV, 
somewhat below the positions of the two calculated peaks. The influence of transfer can 
be seen by comparing the thin curves C12 and Chl3, which do not include the influence of 
transfer, to the corresponding thick curves Ch24 and Ch26, which include the effect. The 
influence is relatively modest but it does shift the calculated peaks towards the peak position 
of the data. 

The S factor for the fusion of ^^C+^^C is illustrated in Fig. [T2j The calculation 'WS 
Ch26' is based on the standard Woods-Saxon potential with radius R^s = 5.352 fm, depth 
Vo = -38.6 MeV, and diffuseness a = 0.63 fm, which was illustrated in Fig. [2f3. It includes all 
26 channels considered in this sections and makes a very good fit to the data above 4 MeV. 
However, a strong hindrance of the data sets in at energies below 4 MeV. The comparison 
shows that the fusion of ^^C+^^C is an example on the fusion hindrance phenomenon which 
has been observed in many medium heavy systems 

The two calculations Ch24 and Ch26 that are based on the M3Y+repulsion potential 
reproduce the data in Fig. [12] at the lowest and also at high energies. The calculation Ch26 
exceeds the data between 4.5 and 5.5 MeV by up to 50%. The calculation Ch24 is in better 
agreement with the data but there are still some deviations. The measured 5* factor is rather 
flat between 3 and 4 MeV and has a modest peak near 4.3 MeV. The calculated S factors 
do not exhibit any sharp peak structures but rise slowly with decreasing energy. 

The two calculations Ch24 and Ch26 represent extreme views of the influence of the 
mutual excitation of the 2"^ and 3~ states in ^^C, i. e., of a quadrupole excitation built on the 
3^ state or an octupole excitation built on the 2+ state. The calculation Ch24 completely 
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ignores it, whereas the calculation Ch26 exaggerates it by assuming that the 2"*" and 3~ 
excitations are independent modes of excitation. In a realistic, microscopic description, the 
mutual excitation of the 2"*" and 3~ states will not be one discrete excitation but will be 
fragmented and spread out over a range of excitation energies. This is indeed the result of 
shell model calculations that are discussed in section VI. 

V. PREDICTIONS OF THE FUSION OF ^'^C+^'^C 

Similar predictions are made for the fusion of ^^C+^^C. The entrance channel potential 
that is used is shown in Fig. |2p. It was obtained with the parameters of the point-proton 
density of ^^C shown in Table I, and the diffuseness associated with the repulsive part of 
the M3Y-I- repulsion interaction was again set equal to the value = 0.31 fm which was 
determined in the analysis of the ^'^C+^'^C fusion data. 

The excitation energies that are considered in the coupled- channels calculations are shown 
in Fig. [131 The solid (red) lines represent the excitations shown in Table I which together 
with the elastic channel give a total of 7 channels. The second 0^ state is Hoyle state. 
The (blue) dashed lines in Fig. [13] are the mutual excitations that are considered in the 
calculations. They include the three mutual excitations: (2+,2+), (2+,3^) and (3^,2+), 
with a one-phonon excitation in each ^^C nucleus, and the two mutual excitations: (2+, 3~), 
with both one-phonon excitations belonging to the same nucleus. There is also a (3^,3~) 
mutual excitation indicated in the figure (with one octupole excitation in each nucleus) 
but it is ignored because the excitation energy is very high. The total number of channels 
consisting of the 7 channels mentioned above (represented by the solid (red) lines in Fig. 
[T3|) and the 5 mutual excitation channels (the (blue) dashed lines in Fig. [T3|) is therefore 
12. The associated coupled-channels calculation is referred to as Chl2. 

The calculation Chl2 includes the mutual excitation of the 2+ and 3~ states, not only 
when the two states belong to different nuclei but also when they belong to the same nucleus. 
Since the description of the (2+,3~) mutual excitation in terms of two independent modes 
of excitation is questionable when the two excitations belong to the same nucleus, it is 
of interest to perform a calculation that does not include that kind of mutual excitation. 
Such a calculation contains 10 channels and is referred to as ChlO. The ChlO calculation 
does contain mutual excitations of the 2+ and 3~ states but only when the two one-phonon 
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excitations belong to different nuclei. The two calculations are compared in Fig. [H] to the 



measured cross sections of Ref. 
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gN31|. The Chl2 calculation exceeds the data substantially 
between 3 and 5.5 MeV, whereas the ChlO calculation is lower and intersects with some of 
the measured peak cross sections. 



The S factors for the fusion of ^^C+^^C predicted by the 
are compared in Fi g. IT5| to the data of five experiments 



ilO and Chl2 calculations 
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four measurements 
measurements 
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32|. The symbols for the 



- |31| are the same as in Fig. [TH whereas the symbol for Spillane's 



32| is indicated in the figure. The calculation Chl2 exceeds all of the data 



points shown in this figure. It is considered an extreme upper limit because it exaggerates, as 
discussed above and in the previous section, the influence of the mutual (2+,3^) excitation 
within the same nucleus. The calculation ChlO, which is much lower, is considered a lower 
limit of the prediction. 

The large difference between the ChlO and Chl2 calculations shown in Fig. [15] is very 
unfortunate. It illustrates the difficulty and uncertainty in predicting the fusion cross section 
due to the poor empirical knowledge of excitations that are built on the and 3~ states. 
We shall eliminate that problem in the next section by applying the results of shell model 
calculations. Another uncertainty, which is similar in magnitude, is due to the empirical 
value = 0.90(7) [22]. This is illustrated in the figure by the thin lines and represents an 
error band of up to ±25%. 

The S factor one obtains in the no-coupling limit is shown by the lowest dashed curve 
in Fig. [151 It is a factor of 2 to 3 below the measured peak cross sections and a factor of 2 
to 5 above the minima of the measured cross sections. In other words, it does not provide 
a background nor does it provide an upper limit of the measured cross sections. It is closer 



to the smooth experimental cross section constructed by Yakovlev et al. 



33|. 



VI. SYSTEMATICS OF THE FUSION OF CARBON ISOTOPES 

The discussion in the previous sections shows that it is possible to develop a fairly com- 
prehensive description of the fusion data for the three carbon systems within the coupled- 
channels approach. A disturbing trend in the analysis of the ^^C+^^C fusion data is the 
need for dividing the data by the factor Sc = 0.843, i. e., for increasing the measured cross 
sections by 19%. However, that is not a very serious problem because the systematic error 
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of the experiment [2j is 15%. 

The need for a renormahzation of the ^^C+^^C data is possibly an experimental problem 
because it is difficult to explain why the measured fusion cross sections for ^^C+^'^C 41 are 
larger than the measured cross sections for ^'^C+^^C Naively, one would expect the cross 
sections for the smaller system to be smaller. This is indeed the ordering that is observed 
in the fusion of carbon isotopes with a thorium target js^. Of course, the problem could be 
in the absolute normalization of ^^C+^'^C data Whatever the explanation is, the basic 
problem is illustrated in Fig. [161 where the measured cross sections for the two systems have 
both been normalized to the Ch26 calculation for the ^^C+^'^C system discussed in sections 
III.B. It is seen that the normalized cross sections are fairly constant but the cross sections 
for the smaller ^^C+^'^C system are larger, on average by the factor 1.02/0.84 = 1.21 (c. f. 
the caption to Fig. [T6|) . The 21% deviation between the two measurements is within the 
large uncertainties of the two experiments, each being of the order of 15 to 30% but it does 
show an unexpected trend. Clearly, it is desirable to have the experimental uncertainties 
reduced in future measurements. 

Another disturbing feature discussed in the previous section is the large uncertainty in 
the cross section predicted for the fusion of ^^C+^^C. Half of the uncertainty concerns the 
influence of the (2'*',3~) mutual excitation which is poorly known experimentally 22|]. One 



35| which predict that E2 



way to overcome this problem is to rely on shell model calculations 

excitations built on the 3~ state and E3 excitations built on the 2'^ state primarily populate 
a 4~ state which is about 4 MeV above the 3~ state. That is similar to the model we have 
used for the mutual (2+, 3~) excitation. However, the B- values of these transitions are only 
about half of the B-values obtained for the corresponding E2 and E3 transitions from the 
ground state to the flrst 2"^ and 3~ states in ^^C. 

The cross sections one obtains for the three fusion reactions discussed in this work, using 
the best flt parameters determined in the previous sections, and the shell model results 
described above, are compared to the data in Fig. [T7]in terms of the S factor. The data for 
^■'C+^'^C have been divided by the scaling factor Sc = 0.843 which results in an excellent 
agreement with the data. The calculations of the fusion of ^^C+^^C and ^^C+^'^C are Chl2 
and Ch26 calculations, respectively, which are similar to those discussed in section IV and 
V. The only difference is that the parameters for the mutual (2"*", 3~) excitation in ^^C have 
been replaced by shell model predictions 35 1. 



20 



The calculated 5" factors for the fusion of ^^C+^^C and ^^C+^^C are essentially predictions 
because the parameters of ion-ion potential and the density of ^'^C were determined in the 
analysis of the ^'^C+^^C fusion data, whereas the density of ^^C was determined by a fit to 
the experimental point-proton distribution. The ^^C+^^C data 0, Q shown in Fig. [TTlform 
a plateau between 2.5 and 4 MeV and has a peak near 4.3 MeV. These features are not quite 
reproduced by the calculation but the fit to the data is not poor; the fit to the data of Ref. 
{4] is stable with respect to the scaling factor Sc and to an overall energy shift AE in the 
center of mass energy, i. e., the P^r data point has a minimum for Sc = 1 and AE = 0. 

Although there is still some uncertainty in the predicted 5* factor for the fusion of 
i2q_|_i2q^ it appears that the calculation shown in Fig. [T7] is consistent with the max- 



ima of the measured peak cross sections. The peak at 4.924 MeV is within the error 



Dand of 



the calculation. The 3 low-energy (triangular shaped) data points are from Ref. 32|]. They 



exceed the calculation tremendously but that is not a problem either because those 3 data 
points are questionable as discussed in Ref. {sgI Is^]- The calculation (with its error band) 
is therefore expected to provide a useful guidance and an upper limit for the extrapolation 
to the low energies that are of interest to astrophysics. 



VII. CONCLUSIONS 



We have explored the fusion of different combinations of carbon isotopes within the 
coupled-channels approach. The fusion cross section is determined by ingoing-wave- 
boundary conditions that are imposed at the minimum of the pocket in the entrance channel 
potential. It turned out that the measured cross sections for the fusion of ^^C+^^C cannot 
be explained accurately at low energies when a standard Woods-Saxon potential is applied 
in the calculations. However, the data can be explained fairly well by applying the shallow 
potential that is produced by the M3Y-|-repulsion potential. This is a characteristic fea- 
ture of the fusion hindrance phenomenon which has been observed in many medium-heavy 
systems. One reason the shallow potential makes it possible to explain the carbon data by 
coupled- channels calculations is that the high-lying excitations become closed channels at 
low center-of-mass energies, not only with respect to reactions but also to fusion. 

The coupled-channels calculations are surprisingly sensitive to mutual excitations and this 
seems to be justified by the comparison to the data. Thus the calculations that include one- 
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phonon excitations alone cannot explain the data so well, neither when the standard Woods- 
Saxon nor the MSY+repulsion potential is apphed. The data are reproduced much better 
when mutual excitations are included and the shallow MSY+repulsion potential is applied. 
The reason is that the direct coupling to these high-lying states is proportional to the second 
derivative of the ion-ion potential and this quantity is much larger for the M3Y-|-repulsion 
than for the standard Woods-Saxon potential because of the strong repulsion. 

The large sensitivity to the excitation of high-lying states is particularly critical when 
the mutual excitations involve the strong octupole excitation of ^^C. The uncertainty in 
the strength of the octupole excitation, and also in the strength of transitions to high-lying 
states that are built on the 2+ and 3^ states, makes it difficult to predict the fusion cross 
sections very accurately when ^^C is involved. However, by adopting the nuclear structure 
properties of the high-lying states predicted by the shell model, one obtains a fairly good 
description of the ^^C-|-^^C fusion data. 

The prediction of the cross section for the fusion of ^^C+^^C, based on the ingoing-wave- 
boundary conditions and the shell model prediction of the coupling to high-lying states, 
exceeds most of the measured cross sections and is consistent with the maxima of the ob- 
served peak cross sections. Thus it appears that the calculation provides an upper limit of 
the measured cross sections. The calculation is therefore expected to provide an upper limit 
for the extrapolation of cross sections into the unexplored territory of very low energies. 

On the experimental side it is very important to get a better absolute normalization of 
the measured fusion cross sections. The current systematic errors are fairly large. This was 
exploited in the analysis of the ■'^^C-|--^^C fusion data by focusing on the shape of the measured 
cross section, whereas the absolute normalization was treated as an adjustable parameter. In 
this way it was possible to calibrate the M3Y-|-repulsion potential and achieve an excellent 
fit to the data. It is of great interest to know whether the renormalization that was applied 
to the data can be justified by more accurate measurements. 
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VIII. APPENDIX: COUPLINGS 



The matrix elements of the couphngs between two heavy ions are generated by matrix 



elements of the static or dynamic deformation amplitudes (see for example Ref. 16|). 
The matrix elements are expressed in terms of the reduced matrix elements (||a;A||)) 

{hM2\ax,\hM,) = (hM.XfilhM^) (14) 

The off-diagonal reduced matrix elements can be extracted from the measured B{EX) values 
according to the formula 

B(£A,/,^« = (5££!M)^K£M>P. (15) 

^ ' ^ ^ 47r ^ 2/2 + 1 ^ ^ 

If the 5-value is expressed in Weisskopf units the relation is 

BiEXU,, = l^(A + 3)FI(/.l|a.l|J.>l^ (j,, 
^ 47r 2/2 + 1 ^ ^ 

The diagonal matrix elements, on the other hand, can be obtained from the measured 
quadrupole moments, 

Qn = )J-^ a/|a2o|//). (17) 

They are related to the intrinsic quadrupole moments Qq and the K quantum number by 

Qn = {IK20\IK) {im\II) Qq. (18) 

For a deformed nucleus with a static deformation I3x the deformation amplitude is axf_i = 
l3\D^Q{e), where e is the orientation of the symmetry axis. In this case the expressions for 
the matrix elements between states \IKM) are 



{hKMlaxolhKM) = ^^^^-±1 {I,M XOlhM) {hK XOlhK). (19) 
This implies (by comparing to Eq. f lT4|) that the reduced matrix elements are 



/2||aA||/i) = /32V2/1 + 1 {hK XO\hK). (20) 

. + 3)]^ 2/1 + 1 
An 2/2 + 1 



B{EXW^. = Mi±M! ^Jl±l i^i^K X^\hK)\ (21) 
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TABLE I: Parameters of the densities that are used in calculating the M3Y+repulsion potential. 
The second column shows the radius R; the diffuseness was kept fixed at a = 0.44 fm and the third 
column is the calculated rms radius. Also shown is the diffuseness a,, that is used in calculating the 
repulsive part of the potential, the value of scaling factor Sc that optimizes the fit to the ^'^C+^^C 
fusion data Q], and the /N. The last 3 lines show the radii that are consistent with the rms 
radii of the point-proton (pp) and matter distributions of ^^C and ^^C. 
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TABLE II: Properties of E2 and E3 transitions in ^^C and ^^C |22[ |. The intrinsic quadrupole 
moments have been extracted from the lowest quadrupole transitions. 



Nucleus 


State 


E^ fMeV) 


Transition 


B(EA) (W.u.) 






"1 







Qo =-19.5 fm^ 


0.570 




2+ 


4.439 


E2: 0^ 2+ 


4.65(26) 


0.570 




0+ 


7.654 


E2: 2+ 0^ 


8.0(11) 


0.236 




3~ 


9.641 


E3: of 3- 


12(2) [22] 


0.90(7) 


13c 


1/2- 







Qo = -17.9 fm2 


0.495 




3/2- 


3.6845 


E2: 1/2- 3/2- 


3.5(8) 


0.495 




5/2- 


7.547 


E2: 1/2- 5/2- 


3.1(2) 


0.465 




5/2+ 


3.8538 


E3: 1/2" 5/2+ 


10(4) 


0.82 




Q is 



FIG. 1: (Color online) Experimental point-proton density of C [19[ is compared to the sym- 
metrized fermi function distribution with diffuseness o = 0.44 fm and radius R = 2.155 fm. 
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FIG. 2: (Color online) The M3Y+repulsion double-folding potentials for different combinations of 
carbon isotopes. They are based on the parameter ar = 0.31 fm and the radius R = 2.28 fm for 
^^C, and i?=2.155 fm for ^^C. The black dashed cures are the pure M3Y potentials. Also shown 
are the WS potentials discussed in the text and the ground state energies of the compound nuclei. 
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5/2" _ 


3/2" 


5/2^ 




3/2" 




1/2" 


1/2" 













FIG. 3: (Color online) Excitation energies of the channels considered in the calculations of the 
fusion of ^^C+^^C. The (red) solid lines are the excitations in projectile and target of the ^^C 
states shown in Table I. The (blue) dashed lines are the 6 mutual excitations that involve the 3/2~ 
and 5/2+ states in projectile and target. 
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23456789 
Ecm. (MeV) 

FIG. 4: (Color online) Enhancement of the fusion of ^^C+^^C relative to the cross section of the 
no-coupling calculation. Shown are the enhancement of the one-phonon calculation (phi), the 
mutual excitation (Chl3) calculation, and the full calculation (Ch26) which also includes one- 
neutron transfer. The calculations are based on a standard Woods-Saxon potential with radius 
Rws = 5.268 fm. Also shown are the ratios of the measured cross sections relative to the 

no-coupling limit. 
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FIG. 5: (Color online) Results of the analysis of the ^^C+^^C fusion data [2] as function of the 
radius of ^'^C. The dashed curves show the /N for the indicated fixed values of ar and minimized 
with respect to the scaling factor Sc- The parameters of the best solution for ar = 0.31 fm are 
shown in Table 1. 
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^^C+^^C, M3Y+rep 



23456789 10 
Ecm. (MeV) 

FIG. 6: (Color online) Ratios of the measured and calculated fusion cross sections for ^^C+^^C. 
The coupled-channels calculations include 26 channels and are based on the MSY+repulsion poten- 
tial using the radius R = 2.17 fm (BF, 5^=1) and i?=2.28 fm (ABF, 5'c=0.843), respectively. 
Also shown is the ratio of Charvet's data [3] and the calculation that uses the radius R = 2.28 fm. 
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FIG. 7: (Color online) Enhancement of the fusion of "'^'^C-I-^'^C [2] relative to the cross sections 
obtained in the no-coupling limit. The curves show the enhancement of calculations that include 
one-phonon (phi) and mutual excitations (Chl3), and the full calculation (Ch26) which also in- 
cludes one-neutron transfer. The enhancement factors for Trentalange's data 4] have been divided 
by the optimum scaling factor, Sc = 0.843. 
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FIG. 9: (Color online) Measured fusion cross sections for ^^C+^^C [sl are compared to the full 
coupled-channels calculation (Ch26, the solid curve). The thin (green) dashed curves are the 
calculated cross sections for a maximum angular momentum of Lmax = 2, 4, 6, 8, 10 and 12. The 
(black) dashed curve shows the no-coupling limit. 
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FIG. 10: (Color online) Excitation energies of the channels considered in the calculations of the 
^■^C+^^C fusion cross section. The (red) solid lines are the states given in Table I. The blue lines 
are mutual excitations considered in the calculations. 
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FIG. 11: (Color online) Enhancement factors of the measured ^^C+^^C fusion cross sections 4jl5| 
relative to the cross sections obtained in the no-coupling limit. The curves show the enhancement 
of calculations that include the 24 (Ch24) and 26 (Ch26) channels described in the text. The 
associated calculations without transfer, Chl2 and Chl3, are shown by the thin curves. 
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FIG. 12: (Color online) 5" factors for the fusion of C+ C [J, l5[ are compared to calculations that 
are based on a Woods-Saxon potential and 26 channels (WS Ch26), and on the MSY+repulsion 
potential with 24 and 26 channels, respectively. The no-coupling limit is also shown. 
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FIG. 13: (Color online) Excitation energies of the channels considered in the calculations of the 
^■^C-I-^'^C fusion cross section. The (red) solid lines are the states given in Table I. The (blue) dashed 
lines are the mutual excitations considered in the calculations. The (3~,3^) mutual excitation of 
projectile and target is indicated but it is ignored in the calculations. 
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FIG. 14: (Color online) Measured fusion cross sections for ^^C+^^C [28|-l3l|] are compared 
ChlO and Chl2 coupled-channels calculations discussed in the text. Th e data of Refs. 
have been shifted in energy {AE = +100 keV 13] and AE = +75 keV 
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FIG. 15: (Color online) Measured S factors for the fusion of ^^C+^^C 
coupled-channels calculations ChlO and Chl2, and to the no-coupling limit. The thin solid curves 
show the error band in ChlO calculations due to the uncertainty in the octupole strength, fS^ = 
0.90 ± 0.07. 
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FIG. 16: (Color online) Ratios of measured fusion cross sections relative to the calculation Cli26 
of the fusion of ^'^C+^^C. The error bars include a 15% systematic error. The average ratio for the 



12 
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C+^^C data |4] is 1.02, whereas the average ratio for the ^^C+^^C data [2] is 0.843. 



> 



o 

CO 



10 



18 



10 



17 



10 



16 



10 



15 



Coupl.-Channels 
No-coupl. 



(data/0.843) 




E,^ (MeV) 



FIG. 17: (Color online) S factors for the fusion of ^^C+i^C [2], ^^C+^^C [4|,l5| and ^^C+^^C 



13/- 



i3r 



2a432| 



are compared to the coupled-channels calculations described in the text. The ^'^C+^^C data have 
been divided by Sc = 0.843 which optimizes the fit to the data. The thin solid curves show the 
error band due to the uncertainty in the octupole strength in ^^C, /S^ = 0.90 it 0.07. The black 
dashed curve is the no-coupling limit for ^^C-|-^^C. 
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